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Abstract. Exact results on particle-densities as well as correlators in two models 
of immobile particles, containing either a single species or else two distinct species, 
are derived. The models evolve following a descent dynamics through pair- 
annihilation where each particle interacts at most once throughout its entire 
history. The resulting large number of stationary states leads to a non- vanishing 
configurational entropy. Our results are established for arbitrary initial conditions 
and are derived via a generating-function method. The single-species model is the 
dual of the ID zero-temperature kinetic Ising model with Kimball-Dcker-Haake 
dynamics. In this way, both infinite and semi-infinite chains and also the Bethe 
lattice can be analysed. The relationship with the random sequential adsorption 
of dimers and weakly tapped granular materials is discussed. 



PACS numbers: 05.20-y, 64.60.Ht, 64.70.qj, 82.53.Mj 

1. Introduction 

Exactly solvable models play an important role in understanding the complex 
behaviour of strongly interacting many-body systems, see [5TJ [45j [8j [57J [60] GH] and 
refs. therein. With very few exceptions, exactly solvable models are largely found 
in one spatial dimension. Exact solutions are important because standard mean-field 
schemes are inadequate in low dimensions where effect of fluctuations is strong, as 
confirmed experimentally. Classic examples are the kinetics of excitons on long chains 
of the polymer chain TMMC = (CH 3 )4N(MnCl3) [41] . and other polymers confined to 
quasi-onc-dimcnsional geometries [501 140] . Systems like carbon nanotubes, which are 
essentially one-dimensional, have been extensively studied; for instance, the relaxation 
of photo-excitations [S3] or the photolumincscence saturation [5S]. What is common 
in all these systems is that particles move diffusively on an effectively one-dimensional 
lattice, and when two particles encounter, one of them disappears with probability 
one. These models can be solved, for example by the well-established method of 
empty intervals [8] 011]. The average particle concentration c(i) ~ f 1 !" 1 [61] found 
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Figure 1. Parameter space of kinetic Ising chains as parametrised by 7 and 5 in 
Qp. The usual Glauber dynamics corresponds to <5 = and the KDH dynamics is 
given by 7 = 25 (dashed line). The zero-temperature line corresponds to 7 = 1 + 5 
(thick line) and the infinite-temperature line to 7 = 0. The full circle at 5 = 1 
and 7 = 2 corresponds to the dynamics analysed here. 



by this method is in complete agreement with the experimental results, unlike the 
mean-field prediction qyrp(i) ~ 

Besides particle systems, one may consider the kinetics of magnetic systems whose 
phenomenology can be captured by kinetic Ising models 29 . We point out that 
such models arc relevant for the description of the slow relaxation dynamics of real 
systems such as the single-chain magnet Mn2(saltmen)2Ni(pao)2(py)2](C104)2 |17j . 
The configurations of the Ising chain are {a} = (o\^n, • • - ol), with the Ising spins 
01 = ±1 attached to each site i of the chain. The master equation describes the 
sequential time-evolution of the probability P({cr};t) of a configuration {a} 

^P({<r};t) = ^[w i (-<7i)P(...-o-i...;t)-c«; i (o-i)P(...o-i...;t)] (1) 

i 

where the . . . means that all other spins arc not affected and the choice of the rates 
Wj((7j) selects a specific model. The most general form of the transition rates for the 
Ising chain which does not impose explicit conservation laws, keeps the global Z2- 
symmetry Cj i-+ — <n and which only takes into account the spin Ui and its two nearest 
neighbours <r»±i, is given by 



1 



7 



where the parameters 7 and 5 are restricted to the intervals < 7 < 1 + 5 and \S\ < 1. 
This parameter space is illustrated in figure [TJ One requires these transition rates to 
obey the detailed balance condition 

Wi(cr i )P oq (. . .(Tj . . .) = W i (-cr i )P eq (. . . - cr, . . .) (3) 

where P cq (. . . a, . . .) = Z^ 1 exp[— "H[cr]/T] is the equilibrium probability distribution 
where T denotes the temperature of the heat bath and H = —JJ2i Vi^i+i IS the Ising 
chain Hamiltonian with an exchange coupling J. Then detailed balance implies 

7 = (1 + 5) tanh(2 J/T) (4) 
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The Glaubcr-Ising model is specified by the choice [55] 

6 = , 7 = tanh(2J/T) (5) 

Here, we shall rather consider an alternative form, proposed by Kimball 38s and Deker 
and Haake [20] (kdh), namely 

= 7 tanh(2J/T) 

2 2 - tanh(2 J/T) { 1 

We remark that along the zero-temperature line 7=1 + 5 the detailed balance 
condition ([3]) is trivially obeyed, since both sides of equation Q vanish (for T = all 
stationary states are absorbing states). 

Physically, the Glaubcr-Ising and kdh models arc very different. This may be 
understood first from the exact relationship between the relaxation time r and the 
spatial equilibrium correlation length £, viz. 

f 2 ; Glauber [23] 
? ' \ 4 ; KDH [201 EH [22 1 J 

which in both models can be established from the exactly calculable time-dependent 
average spin. Second, and more crucially, at T = 0, the Glauber-Ising model has 
only two absorbing states (a feature typical of a relaxing simple magnet), whereas 

the KDH- model has ~ ((1 + y/5)/2j [El [19] absorbing states on a chain of L sites 
(a feature more typical of a relaxing glassy system, such as the Frederikson- Andersen 
model [IB])- While almost any quantity of interest can be computed explicitly in the 
ID Glaubcr-Ising model, see e.g. [29l [51], for the kdh- model only global averages 
[38} [20] or certain global correlators and responses have been found [24] . In the kdh- 
modcl, often the corresponding equations of motion only close for certain classes of 
initial conditions and in the L^oo limit [24] . 

Our work started from an attempt to find more exactly computable averages for 
the KDH-model. For the remainder of this paper, we shall restrict to zero temperature 
T = 0, hence 6 = 1 and 7 = 2 for the KDH-model. Then one can define the dual kink 
variable r\i <G {0, 1} and the rates ([6]) become 

T)i ■= - — , uji(ai) = -(l - cr j (T i _i)(l - <x,cr i+1 ) = 2r)i-ir)i (8) 

Single-spin flips o~i — > —ai are associated to the annihilation or creation of a pair of 
kinks 

rji-x — > 1 - 7?i_i , rji — ► 1 - r\i (9) 

The only allowed transition is thus the annihilation of a pair of kinks. Neither creation 
nor diffusion of kinks are allowed. Therefore, the ID KDH-model at T = is dual to a 
model where the sites are either empty or else occupied by a single immobile particle 
of a single species A, which can undergo pair- annihilation A + A — > + with their 
nearest neighbours and with rate 2. 

Besides being of interest in its own right as an abstract problem in many-body 
physics, there exist several physical motivations for the study of this model: 

(i) Random sequential adsorption (rsa): In the usual USA defined on a regular 
lattice, an atom adsorbed in a site excludes the adsorption of atoms in the 
neighbouring sites. Instead of single atoms one may define RSA of dimers on 
a regular lattice in which the two atoms of a dimer occupy two nearest-neighbour 
sites of the lattice. In one dimension the two versions are equivalent. One 
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then investigates quantities such as the density of particles in the stationary 
state [27]- If wc define the variable rj i that takes the value or 1 according 
to whether the site i is empty or occupied by on atom of the dimer, then the 
relation between this variable and the kink variable rji is simply rj i = 1 — rji 
and the transition rate of the RSA-model is that given by equation (JHJ. A kink 
(r/i = 1) corresponds therefore to an empty site (ry^ = 0) in the RSA- formulation. 
Notice that not any state {rjj} is a possible configuration of dimcrs. However, 
if the initial state is a dimer configuration then any future state will also be a 
dimer configuration. Random sequential adsorption is a commonly encountered 
mechanism when macromolecules or colloidal particles deposit on solid surfaces 
and when gravity effects can be neglected. It goes much beyond the often 
inadequate mean-field-likc descriptions based on the Langmuir equation. Sec 
[541 152] for recent reviews, including experimental tests and applications. 

(ii) Granular matter : the particles making up granular materials rapidly relax into 
one of the very many blocked configurations of these systems. Under the effect of 
a gentle tapping, they may jump from one of these blocked states to another. 
These question has been raised how to formulate a statistical description of 
granular matter in terms of an ensemble of these blocked states. The Edwards 
hypothesis states that all blocked configurations ('valleys') with the same energy 
should be equivalent and that their number can be counted simply in terms of the 
configurational entropy (or 'complexity') [25] . While this proposal seems to work 
very well for mean-field-like models, one of the interest of exactly solvable models 
such as the ones considered here is that because of their large number of stationary 
states, they have a non- vanishing complexity and hence permit a precise test of 
these concepts without appealing to mean-field schemes. Indeed, it has been 
shown that although the Edwards proposal is numerically quite accurate, it leads 
to predictions different from what is found in exactly solved one-dimensional 
systems [H W\ El ECU G23 EH ■ For reviews, see e.g. [3DJHB]. 

Therefore, we shall present here exact derivations for particle-densities and some 
correlators, for two models undergoing a 'descent dynamics' where each individual 
particle only reacts once through the entire history. The first one simply is the pair- 
annihilation model i + i -> 8 + 1 of diffusion-less particles, with the rate ©, and 
dual to the Ising model with KDH-dynamics, see sections 2 and 3. The second model 
describes the pair annihilation /! + £?— i>0 + of two distinct species of diffusion-less 
particles, see section 4. A possible application of the two-species model is related to 
the fact that it could be used to distinguish the two types of kinks in the Ising model, 
namely ^1= A and B, such that one would have the duality correspondences 

(tit— ►m) = (ab — > 00) , (m^m) = (BA^M) 

While the reaction-diffusion processes mentioned above can be treated via the 
empty-interval method, the diffusion-free models we shall be considering here can be 
treated by its converse [27 , which is sometimes called the full-interval method [37 . 
A central quantity will be F n := P( • • .. . • ), the probability of having at least n 

n particles 

consecutive sites occupied by a particle. It turns out that the F n satisfy closed linear 
systems of equations of motion, to be derived from (HjU Here, we shall go beyond 

\ In contrast with the full-interval method as used in |37I|3] , we do not require the condition Fo(t) = 1, 
but shall rather analytically continue the equations of motion to the case n = which will then be 
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the calculation of mere densities. Since in the context of the empty-interval method, 
an efficient way to generate correlators and responses is to consider probabilities of 
pairs of empty sites El HI EH HH H3] , we shall derive correlators by solving the 
equations of motion of the probabilities of pairs of groups of occupied sites, separated 
by a gap. Rather than restricting to rather special initial conditions such that the 
solutions of motion can be guessed by making an appropriate ansatz, wc shall present 
a very general and easy-to-use generating-function method which allows to derive 
systematically all quantities of interest for arbitrary initial conditions. The usually 
considered situation of initially uncorrelated particles is recovered as a special case. 

Binary RSA reactions of immobile point-like particles have also been studied since 
a long time in the continuum. Besides the so-called 'black spheres model', where 
particles closer than a given distance r react instantaneously |42j . other commonly 
studied interactions include 'recombination through tunnelling', where the reaction 
probability w(r) = wq cxp(— r/ro) depends exponentially on the distance r between 
the interacting particles or 'multi-polar interactions' with the form w(r) = Wo(r /rv)~ s . 
In these cases, reactions become possible across arbitrarily large distances which leads 
to a totally different structure of stationary states than the one of the lattice models 
treated in our work. In view of the effective long-range interactions thereby induced 
between all particles, it is not totally surprising that mean-ficld-like approximations, 
such as Kirkwood's 'superposition approximation' [5^1 U (which serves to close the 
infinite hierarchy of equations of motion for densities, pair correlators,. . . and is quite 
similar in spirit to the 'pair-approximation' or '(2,l)-approximation' in the systematic 
analysis in [6]), produce numerically very precise results [HI (55] EH [43] [35] and which 
have been confirmed later by precise analytical studies [HI [10] E] ■ Generically, only 
the case of initially uncorrelated particles (and of equal initial density Pa(0) = Pb(0) 
in the two-species model) is considered, and one obtains for very long times t —> oo 
a very slow decrease of the average particle density pA(t) ~ ln~ d t for the single- 
species model and pA(t) = PB{t) ~ \n~ d ^ 2 t for the two-species model, both with 
recombination through tunnelling in d spatial dimensions (55] [56] [10] and furthermore 
dynamical scaling forms for the correlators [TT|. For multi-polar interactions, one 
obtains pA(t) = Pb{P) ~ t- d /( 2s ~ d ) g3]. However, these models with the Kirkwood 
approximation do not preserve the exponentially growing number of steady states with 
the system size L (of the order ~ 1.618 L for the single-species model and ~ 2.41 L for 
the two-species model) of our discrete lattice models with contact interactions only 
and we do not need to fall back on any approximation scheme. 

This paper is organised as follows. In section 2, we consider the diffusion-free pair 
annihilation model with the simplifying technical assumption of spatial translation- 
invariance. Particle densities and pair correlators are derived exactly, for arbitrary 
initial conditions, via a generating function method. We shall study both the cases of 
an one-dimensional chain [36] [47] [48] [19] and as well the Bethe lattice [26] [44] [2] as 
an example of an effectively infinite-dimensional lattice. In section 3, we generalise to 
the case where spatial translation-invariance is no longer required. In particular, our 
treatment includes the case of an one-dimensional, semi-infinite system so that one can 
study the cross-over from the boundary towards the bulk. In section 4, we illustrate 
further the versatility of our generating function technique by deriving densities [44] 
and correlators of a two-species pair-annihilation model of immobile species A, B of 

used to define Fo(t) by this differential equation. In this respect the full-interval and empty-interval 
methods are quite distinct, although the respective equations of motion look very similar. 
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particles undergoing the reaction A + B — > + 0. Section 5 gives our conclusions. 
The computational details for the single- and two-species models are presented in the 
appendices A, B and C. Appendix D outlines the enumeration of stationary states in 
the two-species model. 
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2. The pair-annihilation model without diffusion 

Consider the diffusion-free pair-annihilation model, consisting of particles of a single 
species A, which sit motionless on the sites i of an infinite linear chain and which may 
undergo the only admissible reaction A + A — > + with rate 2 between nearest- 
neighbour sites. From the master equation (fTJ), with the rates ®, one can derive 
the equations of motion for the averages of interest. In general, any average (/({??})) 
satisfies the evolution equation 

^(/(W))=E/(W)| p (w^) 

= 2 E E /(M) [(! - w-iX 1 - ^) p (- ■ • a - w-i). (! - •••;*) 

-r?i_i7?,P(. r)i...;t) 

= 2 E E [/(• • • ( J - (i - ■ • •) - f({v})]vi-mP(- ■ • w-i, % ■••;*) 
= 2E([/(- ■ ■ (i - m-im (i - *(*)) • • ■) - fm)]m-i(t)m(t)) (io) 

S.i. Density of n-strings 

We shall concentrate on n-strings of at least n consecutive particles • • . . . • (or kinks 

n particles 

in the spin formulation of the kdh- model at T = 0). Their time-dependent density is 
C„(t) :=<»ji (t)i»(t)...»M(t)) (11) 
According to (|lL)p. the C„(i) satisfy the equations of motion 

^Cn{t) = - 2( m (t) m (t) V2 (t) . .. Vn [t)) 

-2(n-l){ Vl (t)r l2 (t)...r ]n (t)) 

-2{r ]1 (t)r l2 (t)...r) n (t)r) n+1 (t)) (12) 
where we used the relations rjf = iji and (1 — r)i)r)i = which follow since rji £ {0, 1}. 

2.1.1. General solution Throughout this section, we shall assume that the initial 
conditions display translation-invariance. Since we shall later follow essentially the 
same approach towards finding the solution, we explain it here, for the most simple 
case, in a little more detail, although the end result for the n-string density C n (t) is 
well-known, at least for uncorrected initial conditions [36 1 f26 | [27 ] l46 | [47 ^ [4H | rT9 ] . 
The evolution equation (|12[) now becomes, for n > 1 

±C n (t) = -4C n+1 (t) - 2(n - l)C n {t) (13) 
at 

To eliminate the last term in (fTU)) . let 

C n (t) = U n (t)e- 2 ^-^ (14) 
so that we are left with 

j t U n (t) = -4U n+1 (t)e~ 2t (15) 
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-it 



1 1 



t = --ln(2s + 1) (16) 



2 2 
so that, for all n > 1 

j-«n(s) = 4u„ + i(s) (17) 
as 

In order to simplify the following computation, and all those shall follow later, we now 
define an auxiliary quantity Uo(s) such that (|T7|) is valid for all n > 0. We shall show 
below that uq(s) does not enter explicitly into any physical observable of interest. 
Equation (|17|l is now solved by introducing the generating function 

F(x,s):=J2^ X " (18) 

n=0 

which in turn satisfies the equation 

£*«.•>-£ - * £ A- - .) d9) 



n=0 n— 1 

which has the general solution F(x, s) = f(x + As). Herein, the last unknown function 
/ is fixed from the initial condition f(x) = F(x, 0). We finally have 

F(x, s) = /(4s + x) = F(x + 4s, 0) (20) 

Since the generating function F is analytic at x — by construction, we can expand 
it according to (fT8|) and find 



F(x, S )=jr / ^±(4s + x) k 



k\ 

V - (4s) fe - 

^ k\ ^nl(k-n)r ' 

k=Q n=0 y 1 



OO OO 



\h— n 



n=0 k=n v ' 

«ro+n(0) 



n. 

n— \m— 



where we also exchanged the order of summation and performed a change of variable 
in the index k. Comparing coefficients of a;, it follows that for all n > 



Un(s)= V-!Wn(0) (4sr (22) 
z — j ml 



m! 

m=0 



Reconverting this to the function U n (t) = u n (s), the n-string density C n (t) is finally 

= E ^T C m +n (0) (e- 2t - 1)™ e^ 1 ' 4 (23) 
ro=0 

since consideration of the i — s- limit shows that u n (0) — C„(0) for all n > 0. 
Clearly, the physically relevant n-string densities C n (t) with n > 1 are independent 
of the initial value Co(0), as it should be. Eq. (|2U)) will become an important initial 
condition when we shall compute correlators below. 
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Figure 2. (a) Time-dependence of the probability C n {t) of at least n 
consecutive occupied sites, for n = [1, 2, 3] from top to bottom and for an 
initially fully occupied lattice with p = 1. (b) Stationary particle-density 
C 1 s t(p) = limt-xjo C„(t) as a function of the initial density p (full line). The 
dashed/dotted lines gives the stationary particle-density on the Bethe lattice for 
several coordination numbers z. 

2.1.2. Specific initial conditions Eq. (|23p is our first result, and gives C n (t) in terms 
of the all initial n-string densities C„(0). In order to illustrate the physical content, 
we consider the case of initially totally uncorrelated particles, with average density p, 
such that, for n > 

C n (0)= P n (24) 

Carrying out the sum in (|23[) is straightforward, hence 

C n (t) = p n e 2p ( e ~ 2t - 1 )- 2in - 1)t (25) 
and shows the characteristic double exponential in the time-dependence [36l [27l l48l 

EBB 02111 

Two values of p have a particular interpretation for the KDH-modcl at T = 0. 

p = h : this corresponds to random initial conditions for the Ising spins so that the 
probability of a kink is ^. 

p = 1: this corresponds to an ordered anti-ferromagnetic initial state <7j = (— l) 1 , and all 
bonds are occupied by a kink. 

In figured, we show the time-evolution of the n=string density. While for n = 1, 
the density rapidly converges towards a finite value in the limit of large times t —> oo, 
larger n-strings (with n > 2) do not survive. In figure [2Jd, we show the dependence 
of the stationary density Ci jS t.(p) = Ci(oo) = pe~~ 2p on the density p. The decrease 
of that density when p — > 1 is intuitively accounted for by observing that for p large, 

§ We point out an important difference with the fairly similar empty-interval method, which considers 
the probability E n (t) of finding n consecutive empty sites, see [8]. If in addition one assumes that 
Eq(£) = 1, the E n (t) obey equations of motion quite similar to (1 1 3 I t and which are consistent with 
Eo(t) = 1 remaining valid for all times 1221 . However, for the diffusion- less case under study in this 
paper, even if one starts from a fully occupied lattice initially, viz. p = 1 and admits Cn(0) = 1 
Vn > 0, one has C (t) = e 2 *^ 1 ^ 2 *) ^ 1. 
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many more particles will find a partner for reaction whereas for smaller values of p, 
more isolated particles will remain. The maximal stationary density is achieved for 
Pm = §, where C 1)Bt .(§) = l/(2e) ~ 0.1839 . . .. 

Usually the initial condition for the RSA of dimers is taken as an empty lattice 
which corresponds to all bonds occupied by a kink. In this case C„(0) = 1 and the 
well-known results for RSA in one dimensions can be off from (|25|) with p — 1, including 
the stationary density Ci(oo) = e~ 2 . It is worth to consider distinct initial conditions 
for the RSA of dimers. For instance, consider a configuration of dimers constructed 
by placing dimers at a chain as follows. Let us denote by AB a dimer so that a 
configuration of dimer will be ABAB$AB®tf)AB . . . where denotes an empty site. 
Walking along the sites of the chain, starting from the origin, if site * — 1 is empty or 
occupied by a B atom the probability of site being occupied by an A atom is p and 
to remain empty is q = 1 — p. If site i — 1 is A, then site i will be B. This defines 
a Markov chain whose stationary solution is such that the density of vacant sites is 
p = q/(2 — q) and the probability of a sequence of n empty sites is C„(0) = q n / (2 — q) 
which we consider to be the initial condition. From equation (|23|) . it follows 



a n 2o 
CJt) = — — exp \2q(e- 2t - 1) - 2(n - l)t] where q = (26) 
2 - q 1 + p 

so that Ci(oo) = qe~ 2q /(2 — q) = pe~ 4:p ^ 1+p ' which is a monotonically increasing 
function of p. The usual initial condition for dimers is recovered for q = 1. 

In order to illustrate what might happen for correlated initial conditions, we 
first observe that the initial values C„(0) = y]JL pi of at least n consecutive initial 
sites are rather cumulative probabilities such that p n are the probabilities of having 
exactly n consecutive occupied sites. As a simple example, consider the case of an 
initial Poisson distribution p„ = A"e~ A jn\ with parameter A, such that 

oo 

c ^) = Hj/ erX ( 27 ) 

l=n 

From cq. (j2"3"]) one immediately has 

rn=0 V V ' / / 

where T(a,z) is an incomplete Gamma function [1]. Although we did not succeed to 
convert this last sum into a known function, in figure^ we show the relaxation for the 
first few values of n and also compare with the uncorrclated case, already discussed 
above in figure The comparison is possible since the parameters p and A can be 
related by Ci(0) = 1 — e~ A = p. The stationary density is shown in figure and 
displays, again, a non-monotonous dependence on the initial condition. This example 
illustrates that because of the large number of stationary states, of the order ~ 1.62 L 
on a lattice of L sites [TH1 Qj|J [53], macroscopic quantities such as mean particle- 
densities depend on the precise form of the initial condition and their correlations. 

A physically motivated example for a correlated initial condition is the two-species 
irreversible pair annihilation A + B — > + of immobile particles A and B. If Ci^it) 
and Ci.s(t) denote the average densities of the particles A and B, it can be shown [44] 
that Ci(t) = 2 {Ci : A{t) + Ci,B(i)) is part of a set of quantities C„(2t) which satisfy 
cq. (|13|) . but with the initial conditions 

C n (0) = l? ( , . B _! ;™P air . , P =^pXpW (29) 
t ^{pA+ Pb)P ; n impair r vr 
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Figure 3. (a) Time-dependence of the probability C n (t) of at least n consecutive 
occupied sites, for n = [1,2,3] from top to bottom and for an initial Poisson 
distribution with parameter A = 1. The grey lines give the corresponding 
results for an uncorrelated initial state with average density p = 0.5. (b) 
Stationary particle-density C 1 s i (A) for the initial Poisson distribution(full line) 
and comparison with the corresponding results for an initially uncorrelated state, 



with density p = 1 



(dashed line). 



and where pa,b are the initial average densities of the A, i?-particles, respectively. We 
shall return to this model in much detail in section 4. 

2.2. Correlation of two n-strings 

We search for exact correlators of a n-string with a m-string, separated by a 
certain distance r, which leads us to compute the averages for observables such as 
r | « • . . . », where the state of the r central sites is unknown. We shall also 

m 

refer to this as a 'string with a hole of r sites'. 

2.2.1. String with an one-site hole In order to extend the methods of the previous 
subsection, we begin with the case of a single-site hole of the form • • .. . » | 1 | > • . . . • 

n rn 

and the average 

Cn, m (t) := (m(t) ■ ■■Vn(t)r ln +2(t) . ..r, n+m+1 (t)) (30) 
For spatially translation-invariant initial conditions, the equation of motion is 

J t C n,m(t) = -2K +1 , m (i) + Cl m+1 {t) + (n + m- 2)Cl m (t) + 2C n+m+1 (t)] (31) 
In appendix A, it is shown that the unique solution explicitly reads 
^LnW = (l-e 4t )C n+m+1 (i) 

oo f<l (a) 
k,l=0 k"'" 
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The correctness of this expression can be directly confirmed by checking the required 
initial conditions and by insertion into pip . The linear system of equations (|3ip is 
known to have an unique solution, see e.g. |34j . 

In the special case of a set of initially uncorrelatcd particles, with average density 
p (or uncorrelated kinks in the KDH- model), we have m (0) = p n+m , C„(0) = p n 
and the (n, m)-correlator for a hole of size r = 1 reads 

Cl m (t) = p n+m [p (e" 4t - 1) + 1] e 2p(e- 2t -l)-2(n+m-2)t (33) 

As expected, a finite value in the stationary limit t — > oo is only found for n = m = 1. 



2.2.2. String with a hole of r sites We now consider the (n, m)-correlator of two 
strings of particles (kinks) separated by a hole consisting of r > 1 sites and define 
the average 

C£,m(*) : = M*) • ' • ^W^+r+lW ■ ■ • Vn+m+r(t)) (34) 

Assuming translation-invariant initial conditions, the equation of motion is, see also 
d 

^c; m (t) = -2[c; +lim (t)+c-; m+1 (t)+(n+m-2)c; M jt)+c;;i im (t)+c;-„i +1 (t)] (35) 

The exact and unique solution, for any given initial condition C£ m (0) is (see 
appendix A) 



Cn,m(*) 



-2(n+m-2)t 

fe,/=0 9=0 cj'=C 

'(e- 2t -l) r (e- 2t 



^n+k+q-q' ,m+l+q' W 



(e 



-2t 



fc! /! ?! 



2 r 



1) 



r+l 



A 



e 2 ^ l C n 



-\-m -\-r 



(*) 



(36) 



(r + l)! 

where we also used the explicit expression (|2"3"|) for the n-string density C n (t). This 
exact expression of the (n, m)- correlator C r n m (t) of two strings, separated by a distance 
r, for any initial condition, is the main result of this section. 

In the special case of initially uncorrelated particles with average density p, one 



has <Z, m {0) 



so that evaluation of the sums in (|3"6")l leads to 



Cn,m(t) ~ C n (t)C m (t) 



-e - 2 («+™- 2 )* (2 P (e" 2t 



l)) r 



-2/ 



1 



■ cxp (2p (e' 2t - 1)) - t Fi (r, r + 1; 2p (l - e" 24 )) 



(37) 



where \F\ is the confluent hyper-geometric function pQ, in agreement with earlier 
results [T5|. For a fixed time i, we observe a factorial decrease with r of the connected 
correlator, which is one of the main differences with what would have been predicted 
for weakly tapped granular matter using Edwards' hypothesis which gives a more 
slow exponential decay Q]|l [SO]- For r> 1 large, we find once more the double 
exponential time-dependence. For the first few values of r, we illustrate in figure [4] 
some aspects of the connected correlator C^ m (t) — C n (t)C m (t) in the stationary limit 
t — > oo. In figure the dependence on p is indicated for some small values of r, 
whereas in figure HJd, the dependence on the distance r is continuously interpolated 
from (|37[) . Clearly, the non-trivial correlations decay very rapidly with increasing r 
such that in the stationary state, only very short-ranged correlations persist. 
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0.030 „ , 1 , 1 , 1 , 1 r 




P 



Figure 4. Stationary connected particle-particle correlator CT i. conn s +' ( a ) 
Dependence on the initial density p, for several values of the particle distance 
r. (b) Continuous interpolation of the dependence on the particle distance r, for 
several values of the density p. The full dots give the values for r integer. 



2.3. Some numerical tests 

In order to check our results further, we also also performed Monte Carlo simulations 
of the pair-annihilation process A + A — > + with immobile particles on a chain 
with L = 65536 sites. An average over 2 • 10 5 histories was performed, in order to 
achieve a statistical precision of the order ~ 10~ 6 . Three different initial conditions 
were considered, namely 

a) ... AAAAAAAAAAA . . . 

b) ... A A A® A A A® A A A® A A A® . . . (38) 

c) ... AAAA$AAAA$AAAA%AAAA$ . . . 

The last two of these represent correlated initial states, with either three or four 
occupied sites followed by a single-site hole. 

From the results of the two previous sub-sections, both the density of the n- 
strings C n (t) as well as the (n, m)-corrclators C r n m (t) can be found. The result of this 
comparison is shown in figure [S^bc, for the three initial conditions defined in ((38)) . 
We also checked, where applicable, that indeed C£ TO (t) = C^ n (4). It is necessary 
to re-scale the 'Monte Carlo time-scale' by 0.553 in order to match the implicit time 
scale in our equations. When this is done, one has from figure [5] a perfect agreement, 
both for densities and correlators, and for both uncorrelated as well as uncorrelated 
initial statesjjj] 

2.4. Densities of n- strings on Bethe lattices 

We now consider a Bethe lattic^H with z links per node, see figure [SJ and we are 
interested in n-string averages C n (t) := (r)i(t)rij(t)r)k(t) . . .) composed of a string of 
n contiguous kinks. The evolution equation for C n (t) involves three kinds of terms 

I Eq. J33H implies that Ch(t) — > as t — > 00 for an initially fully occupied lattice with p = 1. 
1 Following [5], the term 'Bethe lattice' denotes the deep interior of the Cayley tree such that the 
strong surface effects of the latter do not arise for the Bethe lattice. 
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Figure 5. Comparison of the predicted density C n (t) of the n-strings (left panels) 
and of the correlators m (i) (right panels) with Monte Carlo data for the three 
initial conditions (a,b,c) as defined in eq. (08} , of the pair-annihilation model on 
a chain with L = 65536 sites. The lines are the predictions taken from I I23I36H . 
where the time was rescaled by a factor 0.553 throughout. The symbols give the 
Monte Carlo results. 



Figure 6. A small portion of the Bcthc lattice with 2 = 3 neighbours. 



[44l [60] : two terms (z — l)C„+i(t) due to the action of the transition rates on the 
tail or the head of the string, (n — 1) terms C n (t) corresponding to the action of the 
transition rates on two neighbouring kinks of the string and (n — 2) terms associated 
to branching of the string. In the following, we shall make the important assumption 
that these last terms can be approximated by C n +i(t), i.e. that the correlation depends 
only on the number of kinks and not on their relative arrangement. This holds true 
if this condition is realised by the initial conditions. The equation of motion then 
reduces to, for all n > 1 

^n(i) = -2(n-l)(7 n (t) -2(n(z-2) + 2)C n+ i(t) (39) 

We see that the special case z = 2 reduces to the case of the linear chain studied 
above. The unique solution is found by a slight generalisation of the methods used so 
far, as explained in appendix C. The result is Vn > 1 



00 ( ? _ o\k T (n + k + i 
= E ( V 7 Cn +k (0) (e- 2t - if e-^t (40) 



fc=0 



-2 
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This is the main result of this subsection. Because of the asymptotic relation 
x h ~ a T(x + a)/r(x + b) = l + 0(l/x) EQ Eq. (6.1.47)], in the limit z -> 2 one recovers 
the result or the linear chain. 

For the special case of initially uncorrelated particles such that C„(0) = p n , we 
reproduce the well-known form [26j HH [2] 

C n (t) = ^e" 2 ^- 1 )' [1 + (z - 2)p (1 - e" 24 )] ~( 2 +(*- 2 )»)/(- 2 ) (41) 

which in the limit z — > 2 indeed reproduces (|2"5")l . We also see that for z > 3, 
the characteristic double-exponential time-dependence becomes a simple exponential 
behaviour. This means that the double-exponential forms encountered above are a 
peculiar feature of one-dimensional diffusion-free systems. In figure we illustrate 
the dependence of the stationary particle-density C 1 s t(p) on the initial density p, for 
several values of z. Upon increasing z, this becomes increasingly like the mean- held 
result. 

3. Pair-annihilation model: the semi-infinite case 

We now extend the considerations of section 2 to the case when the initial conditions 
are no longer required to be translation-invariant. At the end, we shall use this to 
derive densities and correlators of n-strings for a semi-infinite system. 

3.1. Densities of n- strings with inhomogeneous initial conditions 

When initial conditions are no longer invariant under translation, one has to keep 
track of the first site of the string 

C a - n (t) := (r) a (t)r) a+1 (t) . . . r? a+n _i(t)) (42) 

where n € N as before and a € Z runs over all sites of the infinite chain. The evolution 
equation (TT2"j) reads 

iCo;n(*) - -2(n- l)C a;n (t) - 2C a _ 1;n+1 (t) - 2C a , n+1 (t) (43) 
at 

In appendix A the unique solution of the desired n-string density is derived and reads 

Ca . n (t) = jr Ca ~ k '2!kl ni0) ( e ~ 2t - 1 ) m+fc e- 2(n - i)( (44) 

m,k— 

As seen above in section 2, only for n = 1 a non- vanishing stationary density persists. 
The particular case of homogeneous initial conditions (|23[) is recovered when setting 
C a -fc;m+fc+n(0) = C m+ fc +n (0) and rearranging indices. 

We illustrate the physical content by considering the initial condition 

C a;n (0) = 0(a)p n , 9(a) = j J j ![ ^ ° Q (45) 

which describes a semi-infinite system confined to the positive half-axis a > and with 
uncorrelated particles of average density p. First, our result gives an interpolation 
between the particle density Cb ; i(£) = pexp (p (e~ 2 * — l)) t:: ^f pe~ p right at the 
surface and the form (|25p Coo-,i{t) = pexp (2p(e- 2t - 1)) pe -2 P deep in the 

bulk. Right at the surface, the stationary particle density varies monotonously with p 
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Figure 7. (a) Stationary particle-density C ls ^(p, a) on a semi-infinite chain, 
as a function of the initial density p for the distances a = [0, 2, 4, 3, f ] from the 
boundary at a = 0, from top to bottom (dash-dotted lines for a > 1). (b) 
Continuous interpolation of the dependencies of C 1 s ^ (p, a) on a, for several values 
of p. The full dots give the values of C 1 s t (p, a) for a integer. 



which is no longer the case if one penetrates into the lattice, see figure [7^. From (|44|) . 
we find for the stationary particle density 

a I \k 

°° fe=0 
= pe- 2 ? [1 + (-1) V +a 1 F 1 (a + 1, a + 2; p) /r(o + 2)] (46) 

In figure several aspects of this stationary density are displayed. Figure [7^ 
shows the dependence on p for several distances a from the boundary. The asymptotic 
bulk form is rapidly reached, a few lattice sites away from the surface are sufficient. 
However, the approach towards to bulk is not monotonous. The approach towards the 
bulk becomes considerably slower with increasing initial particle-density. Figure [7Jd 
shows the dependence of the stationary density on the penetration depth a, for a fixed 
value of p. We clearly observe an oscillatory behaviour as a function of a. This is 
natural, since in the stationary state both nearest neighbours of an occupied site must 
be empty, which leads to an effective anti-correlation of the particles. 

3.2. Correlation of two n-strings with inhomogeneous initial conditions 

Analogously to section 2, it is useful to consider first a (n, m)-correlator of two strings 
which are separated by a hole of a single site. Consider the averages 

) (47) 

which satisfies the equation of motion 

jCl.^ m (t) - - 2(n + m - 2)Cl. n>m {t) 

- 2C -l;n+l,m(*) _ 2 Ca;n,m+l (*) _ 4C a;n+m+ i (t) (48) 

which naturally generalises (p?Tj). Similarly, for a hole with r > 2 sites, one has 

C r a . n , m (t) := (V ) (49) 




Figure 8. Stationary connected particle-particle correlator C£ 1 . conn s + from 
eq. H52H . for the initial densities (a) p = 1/2 and (b) p = 1, and for several 
penetration depths a, counted from the boundary. The full symbols give the 
values for r integer. 



and which satisfies the equation of motion, with r > 2 
^C r a . n , m (t) = -2(n + m - 2)C r „{t) 

~ 2Co_i; n +i, m (*) — 2Cq.„ „ l+1 (t) — 2C^. n+1)OT (t) — 2C£. nm+1 (i) (50) 

The solution of the equations (|48l50p is constructed in appendix A, with the result, 
valid for all n, m 6 N, a 6 Z and r > 1 



CJ;n,m(*) = 2 r ( ^ ^ ^ + ^ ^ + ^ ) e 2(r+1)t C a; „ +m+r (t) 



"1<? / \ („-2t -,\k+i+q 



+ e 2(n+m 2)t ^ ^ ^ C^_| ;n+ ^ +g _ g , ro+i . +g , (0) ( ^ j ^ - - - (51) 
£,k=0 <J=0 q'=0 V ' ■ ■ H- 

which is the second central result of this work. 

We illustrate the physical content of this result for initially uncorrclated particles 
of average density p, confined to the right half-plane. Then the initial correlator is 
CJ.„ m (0) = p n+m Q(a) and the initial density was given in P5"j). The stationary 
connected particle-particle correlator is read off from (|5ip . where T(a,z) is an 
incomplete Gamma function [T], 

CJ. lil -C r a; i(oo)C a+r+ i.i(oo) (52) 



fJ 2 c -2p r(l + a,-p) 



r(i + o) 



T(l + r,-2/j) e- 2 T(l + r + a, -p) {-2p) r 



r(l + r) T(l + r + a) T(r + 1 

and depends not only on the distance r between the particles but also on the 
penetration depth a of the leftmost one. This connected correlator therefore describes 
correlations between a particle close to the surface and another one more deeply 
situated in the bulk. In figure |8j the connected stationary particle-particle correlator 
and its factorial spatial decay is displayed for two values of the particle-density. It 
shape strongly evolves with changes of the penetration depth a. 
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4. Two-species annihilation model without diffusion 

We now consider two species of immobile particles, denoted by A and B that can 
occupy the nodes of an infinite chain. The only admissible reaction is the pair- 
annihilation between nearest-neighbour particles: A + B — > + 0, with rate 1. As in 
the single-species model treated above, the number of stationary states on a chain of 
L sites grows exponentially with L and is of the order ~ [y2 + l) , see eqs. (|D.5ID.7|) 
in the appendix D for details. If 77, € {0, A, B} denotes the occupation of the i-th site, 
the transition rates in the master equation (|T|) read 

Ui({r/} -> W}) = ^2 [ S Vi,A$vi+i,B + K,bS, h+u a] <^,0<^ +1 ,0 II £ % ,77< (53) 
4-1. Probability of an alternating n-string ABABAB . . . 

4- 1.1. General solution Since reactions occur whenever AB-pairs meet, we consider 
the probability of having sequences of n particles, consisting of uninterrupted pairs 
ABABAB ... or BAB ABA . . . 

A n := P( ABABAB ... ) = (6 m , A 5 V2tB S m , A 6 V4 , B ■ ■ •) 

n sites 

B n := P( BAB ABA.., ) = (5 nuB 5 m , A 5 V3 , B 6 V4 , A ■ ■ ■) (54) 

n sites 

Starting from the master equation 

^P(M;t) = J2 [P(M;tMM M) P({»?};tM(H W})] (55) 

the evolution equation of A n is 

A A „ = J2 K,aS V2 ,b . . . \p(Wy t t)uH(W} -> M) - PttvY^iM -+ W})} 

= ~ {K,A S V2,B...-5 v[tA 5 Tl ^ B ...)uji({r/}^{ri'})P({ri}]t) 
{vhW} 

= - Yl ( 5 m,A$m,B ■ ■ ■ ~ S^a^b ■ ■ •) [ I| S v 3 ,v' } 

{!)},{!)'} i j<?{i,l + l} 

X (^,A^ +I ,B + ^.B^i+i.A) <^,0<^ +1 ,0P(M;*) (56) 

The two terms in the parenthesis cancel unless j = i or j = i + 1. When it is the case, 
the second term vanishes because of the factor <V,0<V + i,0- ^ remains 

d n+1 

—An = - y^y^^!,A^ 2 ,B ■ ■ • (^,A^ j+ i,B + ^.B^iK+i.^) P ({^}^) 
{„} i=0 
n+1 

= ^^(^ijA^ija.B • ■ • (^?y»,A^77 i+ i,B + ^ .BO~i h+1 ,a) ) 
i=0 

= -B n+1 - (n - l)A„ - A n+1 (57) 

where the first term is due to a transition rate acting on the first site of the string 
(meaning that a B was present in front), the second to transitions in the bulk and 
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the last to an annihilation involving only the last site. Analogously, B n satisfies the 
evolution equation 

d 



dt 



-Br, 



-A n+ i - (n - l)B n - B n+ i 



(58) 



We obtain a set of two coupled differential equations similar to (|13j) . The method 
of solution closely follows the one used in the previous sections and is outlined in 
appendix B. We find for the probabilities of the two n-strings 



A n (t) 
B n (t) 



l Y' fL 

A n (0) - B n (0) + V — (A n+m (0) + B n+m {0)) (e"' - l) 

Z J try,] \ t 



m=0 
oo 



-(n-l)i 



B n (0) - A n (0) + V — (A n+m (0) + B n+m (0)) (e-* - l) 



-(n-l)i 



(59) 



This is the first main result of this section. 



4-1-2. Homogeneous initial conditions We assume that A and B particles are initially 
randomly distributed on the lattice without any correlation among them. If the 
densities are denoted pa and ps (such that 1 — pa — Pb is the initial density of 
vacant sites) then 

"! + (-!)" , 1 -(-!)" [PA 

fp£ 

PA. 



A n ^)=pf +i),2] ^ /2] = {PAP B r 2 

B n {*)=P B n+1),2] pT ] = iP A PBT /2 



2 2 
1 + (-!)" !-(-!)' 



(60) 



where [x] is the largest integer smaller than x. The densities of the n-strings follow 
from (El 



A n (t) = (p A p B )"/2 e -(n-l) t 



l p A +p B \ b 2^papW (e-'-l) + 1 PA — PB 

2 ^JpaPb J 2 ^JpaPb 



-(-1)' 



1 _ 1 PA + Pb \ e -VSw5-(e-'-l) _ 1 PA ~ 

2 ^PaPb ) 2 ^JpaPb 



B n(t) = -{paPb) 



i/2 -(n-l)t 



V 2 y^Ipi / 



+ (-!)« 



2 ^/PaPb ) 2 ^JpaPb 



t Pa - Pb 

2 y/PAPB 

(61) 



and which for n = 1 reproduces the known result [44] . 

In figure [5^, the relaxation of the particle-densities A\(t) and Bi(t) towards 
their stationary values is illustrated. In figure inb, the stationary densities arc shown 
themselves, as a function of the initial parameters pa,b- 



4-2. Correlations between two strings ABAB . . . separated by one site 

We now consider the probability to observe two strings separated by a single-site hole 
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Figure 9. (a) Timc-dcpcndence of the particle-concentrations a(t) = A\(t) 
(full lines), with a = 0.5 and /? = [0.1,0.3,0.5] from top to bottom and 
b(t) = Bi(t) (dashed lines) with a = 0.5 and /3 = [0.1,0.3,0.5] from bottom 
to top. (b) Stationary particle-densities a(oo) = A-^ s ^ (full lines) and 6(oo) = 
B 1 s t (dashed lines) as a function of the initial density pa and for the ratios 
Pb/pa — [0.25,0.50,0.75,1.00] from top to bottom for the A-particles and from 
bottom to top for the B-particles. 



A\ m :=V{ ABABAB . . \ . . ABABAB 

lb, III \ y j I I y 



n sites m sites 



1-/2, b ■ ■ ■ $ri n ,A ( 5))„ + 2,A<5T;„ + 3,B ( 5r;„ + 4, J 4 • • •) ; if 71 impair , . 
(S m ,AS m ,B ■ ■ ■ S Vn ,B S Vn+2iB S nri+3 ,A8 Vn+i ,B ■ • •> ; if n pair 

By definition, the leftmost particle is A (thus the name of this quantity). The first 
block consists of n particles of alternating species A and B. The second block similarly 
consists of m particles of alternating types. The important point is that the last 
particle of the first block is the same as the first of the second one. The particles on 
even (odd) sites are all of the same type. The state of the site in the hole is unknown. 
In the same manner, we also define the quantity 



B! nm := P( BAB ABA...} 1 . . . BAB ABA . 



n sites m sites 

(S Vi ,bS V2 ,a ■ ■ ■ Sti„,b <5i ) „ +2 ,B<V +3 ,A<5 I) „ + 4,-B • ■ •) ; if n impair , . 
(S Vu BS m: A---5 Vn ,A5r ln+2 ,ASr ln+3 ,BSr ln+4 ,A---) ] if n pair 

but here, the leftmost particle is B. By exchanging A and B particles, A\ m is 
transformed into B\ m . The equations of motion are obtained using (|56[) 

^T^n.m = — B n+l,m ~ ( n ~ l)^n,m — 2A n+m+ i — (m — 1)A\ m — Al l m+1 

fo B n,m = - A n+l,m ~ i n ~ l ) B n,m ~ 2B n+m+1 - (m - l)B^ m - B^ m+1 (64) 
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where A n and B n are the correlations of a string of particles without holes previously 
calculated. The solution of these equations is given in appendix B and we find 

^ oo 
j4n + m+l(0) — £? n + m +l(0) + [A n + m+ k + l + -B n + m + fc+l(0)] 



e-< -1) e~ z + l 



\ 



OO 

+ Yl [ A n+k,m+i(0) - B n+kim+l (0)) + [Ai +& , ro+! (o) + Bi +& , ro+! (o)]) 



(■--■) 



and 



Bn im {t) = X 



OO 

(e~* - ^e~* + B„ +m+ l(0) - A n+m+ l(0) + ^ [Ai + m+fc + l + Bn+m + fc+l(0)] 



oo 

+ (V 1 )* [ s i+fe,m+i(0) - Ai +k>m+l {0)] + [A^ +fc , m+; (0) + B n 



+ k,m+l 



(0) 



k,l = 



k + l 



k\l\ 



4-3. Correlations between two strings ABAB . . . separated by r sites 



We now consider the probability to observe two strings separated by a hole of r > 2 
sites and define 



A, 



ABABAB . . , |T~1 , . . ABABAB . 



(67) 



if n impair, r impair 
if n pair, r impair 
if n impair, r pair 
if n pair, r pair 



n sites m sites 

(<5jji,A<5>j 2 ,-B ■ • ■ <5>)„,yt 8r) n+r+1 ,A<)ri, l+r+ 2,B0'r) n+r+ 3,A 
(S, n ,AS, n ,B ■ ■ ■ 0~r ln ,B S, ln + r+1: B5 r , n+r . +2 ,AS rlTl+r . +3 ,B 
^r)2,B ■ ■ ■ S rhl ,A &r) n+r+ i,B°'ri n+T+ 2,A&r) n+r+a ,B 
(S m ,AS m ,B ■ ■ - Sr, n ,B S Vn+r+u ASri n+r+ 2,BSr, n+r+3 ,A 

The leftmost particle is again A (thus the name of this quantity). The first block 
consists in n particles of alternating types A and B. The hole corresponds to r sites. 
The second block consists in m particles of alternating types. The important point is 
that the particle types in the second block is completely determined by the first block. 
If n + r is odd, the first particle of the second block is B while it is A if n + r is even. 
The state of the r sites in the hole in unknown. In the same way, we also define the 
quantity 

■ |"T~| . ..BAB ABA...) 



B^ m := P ( BAB ABA . 



(68) 



n sites 
(S m ,BS V2 ,A ■ ■ ■ § 



\Sr) lt B^f}2,A ■ ■ - SrinjA $7) n+r+1> A$r) n+r+ 2,B&r)n+r+3,A 
(S Vi ,bS V2 .A ■ ■ - 0~r ln ,B 5r )n+r . + 1 ,A0~r l7l+r . + 2,B0~r ln+r . + 3,A 
(S m ,AS, l2 ,B ■ ■ - 0~r ln ,B <V, + r+ i,B<5r ; „ + r+ 2,A<5 J) „ + r+ 3,-B 

In spite of these complex-looking definitions, the equations of motions have the 
relatively simple form, found by using (|56j) 



if n impair, r impair 
if n pair, r impair 
if n impair, r pair 
if n pair, r pair 



-A r 

dt n ' m 

dt ™' m 



-B 



n+l,m 



-A 



n+l,m 



(n - i)A r „ 

(n - l)B r n 



r-1 
n+l,m 

ir-1 
n+1 ,m 



A 



n,m+l 



(m-iK 



n,m+l 



B 



n,m+ 



! - (m - 1)5. 



— /-? r 
n,m n,m+l 



(69) 
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Following the same lines as before (see appendix B for the details), the solution of 
those equation is, for all integers r > 1 



r /{e- t -iy+ 1 (e-'-iy 



(r + 1)! r! 

oo 
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n+m+c+r 

(0) + B 



(70) 
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(71) 



n+m+c+r 

(0) + £?„ +m+c+r 

(0)) 



2 c (e"* - l) c 



c=0 
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which is our last main result. In the stationary state, only the density-density 
correlations A\ 1 and B\ 1 survives. 



4-3.1. Homogeneous initial conditions As before, we assume that A and B particles 
are initially randomly distributed on the lattice without any correlation among them. 
The densities being denoted pa and ps, the initial correlations reads 

'i + (-rr , i-(-i)" rpi 

V PB 



2 2 

i + (-i) m i - (-i) m fi + (-iy 



(72) 



Pb 2 Y Pa 



B r nm (0) = { PA p B ) {n+m)l2 



l + (-l)" , l-(-l) 



2 2 
1 + (-1)™ 1 - (-l) m (1+ (-l)"+ r 



fpl i+r r 

PA 2 

The connected density-density correlations in the stationary state, as given by (f70f 
and (jTTj) in the limit t — » +oo, arc plotted in figure [T0l Qualitatively, the behaviour 
is similar to the one found above in the single-species model. 



5. Conclusions 



We have been studying densities and correlators in interacting-particle models with 
immobile particles and irreversible reactions which have an exponentially large number 
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Figure 10. Stationary connected particle-particle correlator A r ____ „ t (open 

X » X ; COilll ■ . ST/ 

symbols) and 1 . CODn (full symbols) with respect to the distance r between 
the two particles. The different symbols correspond to different values of the 
initial densities. 



of stationary states. These models are relevant in the context of RSA or else for weakly 
tapped granular systems, since they have a non-vanishing configurational entropy (or 
'complexity'). The main results are as follows. In the case of a single type of particle 
with the annihilation reaction A + A — > + 0, we have computed exactly: 

(i) the density (|2"3"|) of the n-strings on the infinite chain 

(ii) the string-string correlator (|36[) of two n-strings on the infinite chain 

(iii) the density (|40j) of the n-strings on the Bethe lattice with arbitrary coordination 
number z > 2 

(iv) the density (|44[) of the n-strings without translation-invariance, including the 
semi-infinite chain 

(v) the string-string correlator (f5Tj) of two n-strings without translation-invariance, 
including the semi-infinite chain 

We did check in all cases that the stated solutions indeed solve the respective equations 
of motion. 

In the case of two species of particles A and B evolving through the reaction 
A + B ->• + 0, we obtained the 

(i) density (|5§f of alternating AB-strings of total length n, 

(ii) the string-string correlator (|70I71|) 

All these results are given for arbitrary physically admissible initial configurations. 
Setting n = 1, they include explicit expressions for the average particle-densities 
and their correlators. In the special situation of initially uncorrelated particles of 
prescribed densities, we reproduce all previously obtained results [36j [26l [23 SH [46j 
S3 S3 [601 OH E] as special cases. 

In general, one observes a rapid relaxation in time, typically described by a 
double exponential form when considering particles on a chain; and by a factorial 
(rather than exponential) decrease of the connected spatial correlators. In the vicinity 
of a boundary, the results deviate strongly from what is found deep in the bulk. 
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In comparison with the KDH model which originally motivated this work, the kink 
densities and correlators under consideration here are distinct from the slowly relaxing 
global observables examined before [201 HI] ■ It would be interesting to see whether 
one might identify a 'dual' to the two-species annihilation model of which global 
variables could be treated in a way analogous to the KDH-modcl. The two-species 
model offers the further possibility to analyse segregation phenomena. 

Acknowledgements: We thank P.L. Kaprisvky for useful correspondence. This 
collaboration has benefited from an agreement between the french COFECUB agency 
and the university of Sao Paulo (agreement Uc Ph 116/09). 
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Appendix A. Calculational details: single-species model 

A.l. Correlator of strings with an one- site hole 

For spatially translation-invariant ini tial c onditions, the equation of motion of the 
averages of the single-site hole ••...• [ 1 | > • . • are given in (|31l) . In the same way 

n rn 

as in eq. (fH| . we set 

C n (t) = u n ( S )e- 2{n - 1]t , < ro (t)=< m ( S )e- 2 ("+— 2 )* (A.1) 
and use again the change of variables (Q~6|) . The equations of motion (|31j) become 

^<m(s) = 2K+i, m (s) + < m+1 (s) + 2(2s + l)u n+m+1 {s)] (A.2) 

and the definition of the u\ m {s) are extended such that they hold true for all 
n, m > 0. We stress that u„q(s) 7^ u n (s) 7^ n (s) are unrelated to the density 
of an unbroken n-string, they are rather a computational device without an obvious 
physical interpretation. Introduce the generating function 

F\x,y,s):= £ -JifU^y™ (A-3) 

n.m—O 

and derive its equation of motion as follows 

d -^-\ x n y m 

— F 1 (x,y,s) = 2 22 [ u n+l,m( s )+ u n,m+l( s )+ 2 ( 2s + 1 ) u n+™+l( s )]^|^J 
n,m— 



(n- l)!m! ^ n ' mx 'n!(m- 1)! 



4(2* +1) £ 



(s) (n + m)! 



,n„ ,m 



(n + m)! n!m! 

n.m— 



-x y 



fe=0 
00 



2 (5 + 1) s) + 4(2 - + ]) £ (I + (A - 4) 



The last term involves the derivative of the generating function F'(x, s) = d x F(x, s), 
see eq. (|18j) . of the n-string density without a hole. It remains 

— F x (x, y, s) = 2 — + — )F 1 (x, y, s) + 4(2s + l)F'(x + y,s). (A.5) 



as \ ax ay 

Equations of this kind are conveniently solved by an appropriate change of variables. 
We choose F 1 (s,x,y) = g 1 {a,f3,s) where a = 2s + x and (3 — 2s + y0 With the 
chosen variables, eq. (|A.5|) becomes an ordinary differential equation 

iL~ 2 ik~ 2 £) pl{x > m s) = ^ (a ' A s) = 4(2s + l)F ' {a + p ~ 4s> s) (A - 6) 

+ To find these, recall that the homogeneous equation (d B — 2d x — 2d y )F = admits any function 
F = f(2s + x, 2s + y) as solution, see e.g. 1341 . 
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where a,/3 only enter as parameters. According to cq. (f2"0)l , F'(x,s) = F'(x + 4s, 0) 
and using eq. (|21l) for F, it follows 

^(a, (3, s) = 4(2s + l)F'(a + (3, 0) = 4(2s + 1) £] ^±^( a + (3) n (A.7) 
os £ — ' n\ 

n=0 

Remarkably, the dependence on s is merely the rather simple explicit one, such that 
the integration readily gives 

g 1 (a, 0, s) = 4(s 2 + s) £ ^±i^ (a + /?)"+ g 1 (a, (3, 0) 

n=0 

= 4(s 2 + .) £ (a + fl- + £ (A.8) 

n— n,m— 

where the integration constant appears as the additional term g (a, (3,0). Next, 
express this generating function in terms of the original variables x,y and s: 

F 1 (x-,y, S )=4( S 2 + S )£^±iM(4 S + x + yr + £ %^ (2s + xT{ 2s + y)™ 
L — ' n\ z — ' n\m\ 

n— n,m— 

k,l ,m— 

oo oo oo l / >-~i \ 

i,l=On=i!m=I v ' v ' 

The identification with the definition of the generating function (|A.3[) leads to 

<M - 4(, 2 + .) £ Mfc+ " + 7 +l(0) (4 S )^ + £ <+* fc y (0) ( 2s )*+' (A.10) 
fc=0 ' *;,;=o 

the unique solution of eq. (|A.2|) . In addition, neither (0), nor Uq m (0) nor uo(0) 
enter into the observables u^ m (s) with n,m > 1, as it should be. After performing 
the changes (|16p and ((HJ), the desired correlator becomes 

Cl m (t) = (e~ 4t - 1) £ "fc+n+m+l(0) 2 fc ( e -2t _ ^ e -2(n + m- 2)i 
fc=0 

+ ^ M n+fc,m+l(°) ^_ 2t _ e -2(n+m-2)i (A.ll) 

Finally, since the first line can be expressed in terms of C n (t), using (|23|) . and since 
C^ m (0) = m (0), the final solution is given by eq. (|32"T) in the main text. 

A. 2. Correlator of string with a hole of r sites 

The equation of motion for the (n, m)-correlator of two strings of particles (kinks) 
separated by a hole consisting of r > 1 sites, for translation-invariant initial conditions, 
are, see also [19] 

jC r n , m {t) = -2[C; i+lim W+C; m+1 (t)+(n+m-2)C; m (t)+^i, r Ji)+C;- m 1 +1 (t)] (A.12) 
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Analogously to previous calculations, set 

C; m W=< m We- 2( " +m - 2)t , (r>l) (A.13) 
and change the time according to (|16[) such that the equations of motion become 

^<,m0) = 2 [<+l,m( s )+<,m+l( s )+< l +l,rn( S )+< l 7m+l( S )] ( A ' 14 ) 

and are extended to all n,m > 0, which defines the auxiliary quantities u^ (s) and 
u 0m( s )- Introduce the generating function 

00 v r (?) 

F r (x,y,s):= J2 ^fV^ m (A.15) 
z — ' nl ml 

n,m— 

whose evolution equation is, for r > 2 



d x n y m 
— F r (x,y,s) = 2 ^ [<+l, m ( s ) + <,m+lO) + <+l,m( s ) +<7m+l( s )] 



n,m— 

_°° „n—l„.m 00 ~,n„.m— 1 



(n- l)!m! ^ ™' mv 'n\(m- 1)! 



1=0 

i = m = l 



= 2 (^ + ^) (^folMH^folM)) (A.16) 

As before for the case r = 1, this is solved by changing variables according to a = 2s+x 
and ft = 2s + y such that the function F r {x,y,s) — g r (a,ft,s) obeys the recursion 
relation 

l- 2 ik- 2 |) ^ ^ = %^ a = 2 {I + 1) ^ * s > ( A - 1? ) 

We express this recursion in integral form for g r (a, ft, s) 

g r (a,ft,s) = 2 J As (JL + JL^g'-ifafra) (A.18) 

In order to evaluate this, we consider first a few special cases. In the special case 
r = 2, the expression (|A.8j) of g 1 {a, ft, s) leads to 

n— 



2s E [ u «+l,m(°)+<m+l( )] ^ 



,i! m! 

n.m— U 
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oo 



d_ 
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; [a + ft) 
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+ g 2 (a,ft,0) (A.19) 



— -— (q + /3) =2^— — ioc + ft) (A.20) 

n— 1 ^ ' 
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n= 1 
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Expanding g 2 (a,f3,0) as 



I ml 



g 2 (a,f3,0) = <m(0) 

z — ' n\ m 

n,m— 

and applying (|A.18|) once more, we find 
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(A.22) 



By further iterating the procedure, we conclude that the general result is 
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(A.23) 



and one can indeed check that this solves the recursion (|A.18[) . Next, we express this 
result in terms of x, y and s: 

F r (x, y, s) = 4'' ( 2^—^ + 4) E ( 4s + x + V)' 
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(A.24) 



in order to identify the coefficients, after an additional rearrangement of the sums 



r+l 



<m(«) ^ + 2 (r + 1) , J W„ +m+ r(s) 



fe+i r-l q 
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(A.25) 



For any fixed initial condition, this gives the unique solution of eq. (|A.14[) . The final 
correlator is obtained after performing the transformations (fTB|) and then (|A.13[) and 
relating the initial coefficients m (0) = u r n m (0) we recover ([35]) in the main text. 
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A. 3. Densities of n- strings with inhomogeneous initial conditions 
The evolution equations for the averages C a - n (t) read 

ic a , n (t) = -2{n-l)C a , n {t) - 2C _ 1;n+1 (t) - 2C a . n+1 (t) (A.26) 
at 

They are solved by letting 

C a;n (t) = u a . n (s) (A.27) 

and change time according to (|1 6|) such that the equations of motion read, extended 
to all n > which defines u a -o(s) 

7-« a; n(s) = 2u tt _i ;Tl+ i(s) + 2it a:n+ i(s) (A. 28) 

ds 

Assuming for the time being a spatially infinite system, introduce the generating 
function 



oo oo 



F(p, X ,s):=Y^ E (A.29) 



n< 

n— a— — oo 



The evolution equation becomes 

OOOO f \ OOOO / \ 

a— — oo n= a — — oo n— 



OOOO / \ OOOO /\ 

n— a— — oo n— 1 a— — oo 

= 2(1 +p)—F(p,x,s) (A.30) 

where the variable p appears only as a parameter. Consequently, and taking the initial 
condition into account, the generating function has the form 

F(p, x, s) = F(p, 2(1 + p)s + x, 0) (A.31) 

Expanding according to (|A.29|) and using Newton's multinomial formula 



F( P ,x,s)= e E^^° (2(1 

a— — oo n— 

oo oo /„\ oo I 
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By comparison with eq. (|A.29[) , we find 



.,„(*)= £ Ua - k ' m + h + n ^ (2s) m+k (A.33) 



to! kl 

,k=0 



and finally obtain (|44|) in the main text. 

The particular case of homogeneous initial conditions (|23[) is recovered when 
setting C a _fc ;r „ +fe+ „(0) = C m+fe+rl (0) and rearranging indices. 
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A. 4-- Correlation of two n-strings with inhomogeneous initial conditions 

A.4--1- String with an one-site hole Analogously to sections 2 and A. 2, it is useful 
to consider first a (n, m)-correlator of two strings which are separated by a hole of a 
single site. Consider the averages C\. n m (t) which satisfy the equation of motion 

^Cl, n , n {t) = ~ 2(n + m - 2)C l a . nim {t) 

- 2 Co-l;n+l,m(*) ~ 2Cl. nm+1 {t) _ 4C a;n+m+ i (t) (A. 34) 

natural generalisation of (|3"Tj) . Define C\. nm {t) = u^. n m (s)e~ 2 ^ n+m ~ 2 > t and also use 
(1161) to obtain 



ds 



4n,m( s ) = «o-l;n+l,m( s ) + W o;n,m+l( S ) + 2 ( 2s + l)«a;n+m+l(s) 



(A.35) 



extended to all n, m > 0. The next step is the definition of the appropriate generating 
function 

oo oo 1 / \ 

F\p, X ,y,s):= J2 E ^^-p a x n y m (A.36) 

a— — oo n,m— 

which can be shown, analogously to what has been done above, to satisfy the equation 
ldF 1 (p,x,y,s) _ ^dF 1 (p,x,y,s) | dF 1 {p,x,y,s) 
2 



P- 



■+2(2s + l)F'(p,x + y,s)(A.37) 



ds dx dy 

where F = F(p,x,s) is defined in (|A.29|) and F' = d x F. To solve this, we perform 
the change of variables 

G 1 (p,a,/3,s):=F 1 (p,x,y,s) , a := 2(1 + p)s + x + y , /3 := x - py (A.38) 
and find the simplified equation 
l_dG l {p,a,P,s) 
2 



= 2(2s + l)F'(p, a - 2(1 + p)«, s) = 2(2s + l)F'(p, a, 0) (A.39) 



where in the last step eq. (|A.31j) was used. Integration gives 

G 1 (p, a, (3, s) = 4s(s + 1)F' (p, a, 0) + G 1 (p, a, /3, 0) 
and from the initial condition and the definition of G 1 , we also have 

\ 1+p 1+p 

Inserting and converting back to the original variables, a straightforward but just a 
little tedious computation leads to 
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and from which one reads off the coefficients 

(2sY+ k 1 
t\ k\ 



u a;n,m{ s ) — /j| r,| u a-l ;n+£ ,m+k (0) 



e,k=o 



2s ( S + 1 )EtE( / ) «-««+k+m+l(0) 



Aj=0 

= ^Tr u U+Am+fc( ) +2s(s + l)u o; „ +m+1 (s) (A.43) 

£,fc=0 

At this point, it is an useful exercise to check that this is indeed the unique solution 
of the equation of motion (|A.35p . 



A. 4-2. String with a hole of r sites The general (n, m)-correlator of two strings 
C T a . n m , separated by a hole of r sites, satisfies the equation of motion, with r > 2 

^C r a . n , m (t) = -2(n + m - 2)C r „{t) 

— 2C^_ 1 . n+lm (t) — 2C£. n)Tl+1 (t) — 2C a . n+lm (t) — 2C a . nm+1 (t) (A. 44) 

This may be simplified as usual by defining C^. n m {t) = u r a . n m {s) e - 2 (™+ m - 2 )* an d 
using (JTBJ), we obtain, for r > 2 

~^ u a;n,m( s ) = u o-l;m+l,m( s ) + U a;n,m+l( S ) + U a;n+l,m( S ) + U a;n,m+l( S ) (A. 45) 

In principle, we could again define a generating function and solve the corresponding 
differential equation. However, the great similarity between the results of the previous 
sub-section of this appendix allows us to write down immediately the expected solution 

u a;n,mW 0\ ]A L < L / *-* :n+t+a-a' ,rn+k+a' W 1 „' I | 

«,fe=0 ' ' q=0q'=0 \ / 

s r+l \ 

" + 2 ( r+1 y ) Ua:n+m+r(s) (A.46) 

valid for n, m € N, a £ Z and all integers r > 1. Here, we have also used the solution 
(|A.33j) for string-densities without a hole. From this, one readily recovers (|51[) in the 
main text. 

Indeed, it is a straightforward matter to verify that the solution (|A.46|) indeed 
solves eq. (|A.45[) and it obviously satisfies the required initial condition for s = 0. 
Hence we have found, for any fixed initial condition, the unique solution of (|A.45[) . For 
a rapid check, note that in the case of spatially translation-invariant initial conditions, 
we recover our previous result (|A.25|> . 



Appendix B. Calculational details: two-species model 

B.l. Densities of n-strings 

The two string densities A n and B n obey the equations (|57|58[) , which are quite similar 
to (fT3| . The method of solution close follows the one used in appendix A. Perform 
the change of variables 

A n (t) = u n (s)e-^-^ , B n (t) = v n (s)e-^-^ ; s := e~* - 1 (B.l) 
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so that the equations of motion become, for all n > 

7-«n(s) = u n+ i(s) + v n +i(s) , -r-Vn(s) = v n+ i(s) + u„+i(s) (B.2) 
as as 

These two equations can be decoupled by the further change of variables 

X n (s) := u n (s) - v n (s) , Y n (s) := u„(s) + v n (s) (B.3) 
for which the evolution equations read 

^-X n (s) = , ^-Y n (s) = 2Y n+1 (s) (B.4) 
ds as 

The first equation gives trivially X n (s) = X„(0) while the second one reduces to (|T7|) . 
up to a time-rescaling so that we have Y n (s) = C n (s/2), already analysed in section 2, 
with the proper identification of the initial condition. Adapting (|23p , one has 

-i 1 1 00 9 m 

u„(s) = -(X n (s) + Y n (s)) = -X n (0) + ^ Y m+n(0) ( e - ( - l) m 

i 1 1 °° 9 m 

»»(*) = 2 ( y "( fi ) ~ X "( s )) = -2 X "(°) + 2 H ^T y ™+«(°) - l T ( B - 5 ) 

7Ti=0 

and this finally leads to the result stated in eq. ([59]) . 

As a preparation for the computation of the correlators below, define the 
generating functions 

00 n 00 n 

F(x,s) := Y]un(s)—r , G(x,s) := Y%„(s) — r i 

ra=0 n=0 
oo n oo n 

X(x,s):= ^X„(s)^- , Y(x,s) := 5>„( S A (B.6) 

n=0 n=0 

which are to be found by solving the differential equations 

dF_ _ dF_ + dG _ dG dX_ _ Q dY _ ^dY_ 
ds dx dx ds ' ds ds dx 

Hence the equations for X(x, s) and Y(x, s) decouple and we recognise those already 

treated in section 2 and appendix A. We are thus back to (|59]) . 

B.2. Correlations between two strings ABAB . . . separated by one site 
The equations for two n-strings separated by hole of r = 1 site read 

~^A nm = — B n+ i^ m — (n — l)A nm — 2A n + m +i — (m — l)A nm — A nm+1 

fo B n,m = - A n+l,m ~ ( n ~ ^) B n,m ~ 2 - B n+m+l ~ (m - l)B^ m - Bj l m+1 (B.8) 

where A n and B n are the correlations of a string of particles without holes previously 
calculated. Perform the habitual change of variables s = e~* — 1 and 

< m « = < m (*)e- (n+m - 2)t , < m (t) = < m (a)e-("+™- 2 )* (B.9) 
while for the non- interrupted A n and B n , use again (jB.lj) . This leads to 

^"n.rnO) = «n+l,m( s ) + 2 ( s + l)u n+m+1 (s) + U^ m+1 (s) 

T- V l, m (s) = <4 +1 , m (s) + 2(s + l)« n+m+ i(s) + v^ m+1 (s) (B.10) 
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This set of equations cannot be decoupled at this point by defining linear 
combinations of them as in (|B.3|) . We therefore introduce the generating functions 



(B.ll) 



1 1 m I 



F 1 (x,y,s):= V u\ m {s)— — , G 1 {x,y,s):= V i£ (a)- 



n,m— n,m— 

and recall the definition (|B.6[) of s) and G(x, s) for the correlation without hole 
Multiplying the evolution equations by x ^ \ and using the fact that 

Un+m+l — j — : 



L»*+i^T- = ^ Uk (fc-1)! 



n,m=0 fc=0 fe=l 

we obtain (with F'(x, s) = d x F(x, s) and G'(x, s) = d x G(x 1 s)) 
OF 1 dG 1 n/ . 3F 1 



E 

.fe=0 



k\ 



(B.12) 



9s <9x 
dG 1 dF 1 



+ 2{s + l)G'(x + y) + 



dy 
dG 1 



(B.13) 



ds dx dy 
Letting now 

X r {x,y,s) = F 1 (x,y,s) - G 1 (x,y,s), 

Y\x, y, s) = F\x, y, s) + G 1 (x, y, s) (B.14) 
the set of equations decouples 

dsX 1 = {-d x + d y )X 1 + 2(a + l)X'(x + y, s) 

dsY 1 = (d x + d y )Y 1 + 2{s + l)Y'(x + y, s) (B.15) 

where X(x,s) = F(x,s) — G(x,s) and Y(x, s) — F(x,s) + G(x,s). The solutions of 
the homogeneous PDEs are 

X\x,y,s) = X\s-x,s + y) , Y^x, y,s) = Y\s + x, s + y)(B.16) 

In the following, denote a± := s ± x and j3 := s + y and perform the following change 
of variables 

X r {x,y, S )=x\a.,p,s) , Y x (x, y, s) — y x (a+,fi, s) (B.17) 

Since 

/ da_ da- 



d s x 1 -(-d x +d y )x 1 = dsx 1 - 

and 



\ ds dx 



t) da - xl+ {^ ~ %) df3xl = ^t 5 - 18 ) 



X'(x + y,s) = —J2x n (0) 



(x + y y 



n=0 

/ ■ \n-l 



(x + y) r 



n=l 

oo 



71=0 

oo 



(n-1)! 
(i + ?j) n 



n! 

(-a_ + 0)* 



oo 



n oro 

^ — ' nl to. 



(B.19) 



n,rn— 
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because X n does not depend on s according to (|B.4|) . the equation of motion for x 1 
becomes 

a n _f3 m 



d s x 1 (a^,f3,s) = 2(s+l) J2 (~l) n X n+m+ i(0)- 
This is readily integrated and gives 



n,m— 



x 1 {a-,p,s) = (s 2 +2s) V (-l)™X„ +m+1 (0)^--+x 1 (a_,/3,0) 

^ — ^ rj.l 777.' 



n,m— 

The last term is written as a series 



X 1 (a.,P,0)= (-l)"^n,m(0) 



n\ ml 



(B.20) 



(B.21) 



(B.22) 



The presence of the oscillating factor (—1)™ is required for X\ TO (0) to be equal to the 
initial value of X 1 , as can be verified a-postcriori. Its origin can be found in the fact 
that a_ = s — x leading to a factor (— x) n in the series when t and thus s vanishes. 
The generating function X 1 is now 

X 1 (x,y, s) = x 1 (s - x,s + y,s) (B.23) 

(s-x) n (s + y) m 



(-1)" [(s 2 + 2 S )X n+ro+1 (0) + Xl m (Q)] 



J.I 777 I 



71,777 — 



Y, (-l)" +fe + 2s)X n+m+k+l+1 (0) + Xl +k<m+l {Q)] n , l fc ,J,^ 

n.in.k.l— 



E 



7i,m,fc,2=0 



(,s 2 + 2 fl )X B+m+fc+I+ i(0)L£j- + (-l) fc Xi +fc>m+i (0) 



fc!Z! 



kill 



x n y m 

nl ml 



E 

n.m— 



( S 2 + 2 S )X„ +m+1 (0) + £(-1)*^ 



+ fc,m+Z 



(0) 



fc,i=0 



„fc+2 



fc!Z! 



77 I 777 I 



where we have used the fact that 



E 

k,l=0 



( -s) k s l 

kill 



n=0 fe=0 V ' n=0 



(B.24) 



Similarly, the second equation (|B.15|) gives 

9 s y 1 =2(.s + l)£r„ +1 ( S )^±^ 



oo 



2( S + 1) Y Yn 



n,i=0 



2( S + 1) J] y « 

n,Z=0 
oo 

2(a + 1) £ Y n +m+k+l+l (0) 



(2 3 ) ; (x+ y r 

n!Z! 

(2s) l (a+ +(3- 2s) 1 



n,777.,fc,i=0 



n!Z! 

(2s)'a r |/3" l (-2s) fc 
n! m! fc! /! 



= 2( s + i) £ y, 



71+771 + 1 



(0) 



7 1 777 ! 



(B.25) 



71,771=0 
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where we have used (|B.5|) . Again the s-dependence disappears and the integration 
leaves 



V 



The last term is written as a series 



\a+,l3,s) = (s 2 + 2s) J2 *Vm+i(0)- 



n\ ml 



y 1 (a+,0,O)= CP 1 )^! 

n.m— 



n! m! 



V(a + ,/3,0) (B.26) 



(B.27) 



The generating function Y 1 is thus 
^foj/is) =y 1 (s + a; ! s + j/,s) 

OO 

= £ [( s 2 + 2 S )y; l+m+1 (o) + r n 1 im (o)] 

n.m— 



(B.28) 



(s + x)"(s + y)' 



J.I 7TJ.I 



£ [(s 2 + 2s)F n+m+fe+m (0) + Yl 



s k+l x n y m 



n.m,k,l—0 



E 



k.l=0 



s k+l 
kill 



x'"y"" 



Taking the sum and the difference of (|B.24|) and (|B.29j) . we recover the generating 
functions F 1 and G 1 and then the coefficients u\ l and v\ l by identifying the 
coefficients of x , y , : 

nl ml 

1 S 2 + 2s 

\, m («) = 2(^™.™( S ) + y ",m( s )) = 2 



-Xn+m+l(0) + V, 



(0) 



OO 

o X/ (( _1 ) fe ^ri+fc,m+i(0) + ^n+fe,m+i(0)) 



and 

1 s 2 + 2s 

W «,m( S ) = g ( y ™,™( S ) - X n, m {s)) = ~ 



kill 
(B.29) 

(2s)* 



—X n+m+ i(Q) + Y n + m j,-k+l (0) 
fc=0 

^ ((-i) fc+1 ^ + ^ m+i (o) + r„ 1 +fe , m+/ (o)) 



k.l=0 



kin 

(B.30) 



Applying now the transformations s = e 1 — 1 and (|B.9[) . we finally obtain 



Ai,m(*) 



/ 00 2 fe (e~ t -l N ) 
(e-* - 1) (e-* + 1) ( A„ +m+1 (0) + ^ F„ +m+fe+1 (0) — i— ^ i 



fc=0 



£ ((-i) fc Ai +fe . m+i (o) + y„ 1 +fe , m+z (o)) 



fci=0 



fc!2! 



-(n+m— 2)t 



(B.31) 
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and 



B n,m(*) = 2 



f - i)(c-* + 1) -x„ +TO+1 (o) + ^y„ +ro+fe+1 (o) 



2 fe fe-*-l) 



fc! 



+ (( _1 )' C+1 ^n+fc,m+i(0) + ^+fc,m+i(0)) 

fe,i=0 

as stated in the text. 



fc!/! 



3 -(n+ra-2)t 



(B.32) 



5. 5. Correlations between two strings ABAB . . . separated by r sites 

The equations of motion for the probability to observe two strings separated by a hole 
of r > 2 sites are 
d_ 

dt 
d 



A r = — B r — in — \\A r — 4 — A r 

n,m n+l,m V / n,m n+l,m n,m+ 



1 ( m — l)Ai,m — ^n, m +l 



.4 



n+l,m 



(n - 1)-BI ro - 



! - (m - 



— R r 

n,m+l 

(B.33) 



A first simplification is obtained by making the change of variables 

>W*) = «U(«)e- (n+m - a)< ' = < ro (s)e- (n+m - 2) * ; » 

and we have 



1 (B.34) 



^ s u n,mi s ) — v n+l,m( s ) + U n+l,m( S ) + U n,m+l( S ) + u n,m+li s ) 

-J^ V n,m( s ) = <+l,m(s) + <+l,m( s ) + *Cn+l( s ) + <,m+lO) ( B - 35 ) 

As before, we assume these equations to be valid for all n, m > and emphasise that 
u n o 7^ U n 7^ u o n an d similarly for the v. We then can introduce the generating 
functions 



F r (x,y,s) 



3C 



7/ 



s) 



x n y m 

77. ' 777. ' 



G r (a;,y,s) = ^ < m (s) 



n!m! 



(B.36) 



n,m— n,m— 

Multiplying the evolution equations by , we obtain 

d s F r = d x G r + d y F r + (d x + d y )F r ~ 1 
d s G r = d x F r + d y G r + (d x + d y )G r ~ l 

Letting now 

X r (x, y, s) := F r (x, y, s) - G r {x, y, s), 
Y r (x,y,s) :=F r (x,y,s)+G r (x,y,s) 
the set of equations decouples 

d S X r = (-d X + dy)X r + (d X + 8y)X r - 1 
d S Y T = (d X + 8y)Y r + (d X + dy)Y r ~ 1 

As before in the case r = 1, and following standard techniques [34], we denote 
a± := s ± x and f3 := s + y and we perform the following change of variables 

X r (x,y,s)=x r (a.,0,s) , Y r (x,y, s) = y r (a+, p, s) (B.40) 



(B.37) 



(B.38) 



(B.39) 
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Since the relation (|B.18j) holds also for X r , the PDE for x r is 

d a x r {a-,p,a) = (d x + d y )X r - 1 = (-d a _ + d^- 1 (B.41) 
The case r = 2 is obtained from (|B.21[) : 

9 s a; 2 (a_,/3, s) 

OO j-^ r) *yy\ 

= {-d a _+d l3 ) ]T [( S 2 + 2 S )(-irX n+m+1 (0) + (-ir^. m (0)]^ T (B.42) 

n,m— 

= [ 2 ( s2 + 2s)(-l)"X„ +m+2 (0) - (-ir +1 ^ +lim (0) + (-l)«X^ m+1 (0)] 

n,m— 

which leads after integration to 

x 2 (a_,/?,s) = x 2 (a_,/3,0) 



n\ ml 



(B.43) 



r / 3 \ 

+ 2f^+,s 2 jx (1+m+2 (0) + Xi +l!m (0) S + ^ m+1 (0).s 

n,m— L \ / 

where 

£ 2 (a_,/3,0) = ]T (-l)"^ ro (0)^ (BAA) 

n,m— 

We can infer the solution for a general r : 



.t>_,/3,,s) = E 2 M r— w + - (-i)"^«+«+r(o)- . 

^-^ V (r + 1) r! / n! m! 

+ E <-i)"EE ? p 1 r S !^;U i -p(<')Tr^r om»j 

n,m=0 i=0 p=0 ^ v 

Expressing now this quantity in terms of the original variable x and y leads to 
X r (x, y, s) = x r (s ~x,s + y,s) 

i (s - x) n {s + y) m 
' rrt.l 



E ^ 7— W + T (-l) n ^n+m +r (0) 

' \ fr + 1)! r! / n!m. 

n.rn— v 7 

V r ii"VV - .v r -' mi s ' (*- a n (* + y) r 

n,m=0 i=0 p=0 ^ v ^' 

/ s r + 1 s r \ s aJrb x n y m 

/ , 2 ' 7 TT7T + ~7 (- 1 ) a ^n+m+a+6+r(0)— j j rT 

\(r + l)! r! / n!m!ao! 

a,fa = 

r_1 | /I J „a+b„n..m 

Ef if W x r ~' mi y 

OO 



E 2r ( (^+1)! + H ) *™+"^°) 



r-l 



7J| g i g a+b x nym 

E pji^p+a.nvf I-p+6(°)][ „, o! m | 6 , ( R46 ) 



Z=0 p=0 
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where we have used (|B.24[) . 

The same procedure has now to be applied to the calculation of y r . Similarly to 
(JHUHJ), the PDE for y r is 

d s y r (a+,p, s) = (d x + d y )Y r - 1 = (d a+ + fy)^" 1 (B.47) 

The case r = 2 is obtained from (|B.26[) : 

d s y 2 (a + , (3, s) = (d a+ + d (j ) [(* 2 + 2s)Y n+m+l {Q) + ^(0)] ° 



.i\ ml 

n,m—0 



oo 



n am 

= £ [2( s 2 + 2, s )r„ +m+2 (o) + y^ 1 , r ^o) + y„V l+ i(°)]^^- 48 ) 

n,m— 

which leads after integration to 

y 2 (a+,/3, S ) = y 2 (a + ,f3,0) 



n,m— 

where 



oo r / 3 \ 



rrt.l 



(B.49) 



y 2 (a + ,M)= £ l£ TO (0)^L. (B.50) 



n,m— 

We can infer the solution for a general r : 



n.m— 



F> + ,/3, S )= £ 2 ^( ___ + _) K„ +m+ ,(0)- 



7 1 ml 



+ E EE !( , L ^iC+i-pt )^ (B.51) 

n,m=0 1=0 p=0 * v t ' 

Expressing now this quantity in terms of the original variables x and y leads to 
Y r (x, y, s) = y r (s + x,s + y, s) 

s r+1 s k \ v fn ^(s + x) n {s + y) m 

11 ' m.\ 



n,m— 

oo r— 1 / 



V w 11 Y k ~ l for ( s - x ) n ( s + vY 

Z^Z^ „! (/ _ v )\ n+p,m+l-p\ I p „l m l 
n,m=0 (=0 p=0 ^ v ^' 



n,m,c— 



E 2r ( (7TT)! + H ) W^(o) 



(2s) c x™^ 



+ E EEp ! (;_p)i y n+P+a,m+;- P + & (°)7[ , , , fc , (B.52) 

Finally the original generating functions F r = \{X r +Y r ) and G r = \{-X r +Y r ) 
can be formed and the coefficients u r n m and v r n m extracted from their series 
expansions. We find 
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(r + 1)! r! 



^n+m+r(0) + Y 7 



n+m+c+r 



(0) 



(2*) c 



(B.53) 



oo r— 1 / 



a,b=0 i=0 p=0 
1 

OO 7*— 1 I 



(r + 1) 



c=0 



n+m+e+r (0) • 

c! 



(2s) 



J+a+b 



~*~ ZZ S pi f 7 _ v )\ ^ + -^n+p+a,m+i-p+b(0) + ^n+p+a,m+7-p+b 7! a ! frl 

a,6=(| 1=0 p=0 F ' ' Ph 



Re-expressing this in terms of the desired correlators, we arrive at eqs. (|7QI7ip in the 
main text. 



Appendix C. Calculational details: single-species model on the Bethe 
lattice 

The n-string averages C n (t) := (r]i(t)r)j (t)r)k(t) . . .) satisfy the equations of motion 
d 



df 



C„(t) = -2(n - 1)C„ - 2 (n(z - 2) + 2) C n+1 



(C.l) 



To solve these, let 

C n {t) = u n {s)e-^ n -^ (C.2) 
and with the change of variables (fT6| we have the equations of motion, for all n > 
d 



ds 



i„(s) =2(n(z-2) + 2)u„+i(a) 



As before, introduce the generating function 

U n (s) „ 



F(x,s) :=£ 



(C.3) 
(C.4) 



n=0 



Applying the equation of motion (|C.3j) gives the equation 



d 



gF(x, S )=2j2Hz-2) + 2] 



u n +i(s) 



2 £[(n-l)(*-2) + 2] 



n=l 



(n-l)f 



(* - 2) £ «n(«) (^yy + 2 £ «„(*) 



n=l 



(n-1)! 



F(i,a) 



(C.5) 



This PDE is solved by separation ansatz F(x, s) = f(x)g(s) so that eq. (|(J.5|) becomes 

g'(s) 



9(s) 



{z-2)x—— + 2—— 
fix) fix) 



(C.6) 
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leading to g(s) 



-As 



and an equation for f{x) 



2(z-2)xf"{x)+if'(x) + Xf(x)=0 
that can be recast as a Bessel equation whose solution is 



/o 



2\x 
z-2 



(C.7) 



(C.8) 



where J and Y are the Bessel functions of the first and second kinds. The Bessel 
function z~ v Y v {z) is always non-analytic at the origin z = 0. Such a behaviour 
is incompatible with the analyticity built into the generating function F(x 1 s) by 
construction. Hence we must have C\ = for any A. Moreover, negative values of A 
are not allowed since they would induce a complex generating function. It remains 



F(x,s) 



dXBxx^'^J ,-4 




-As 



(C.9) 



Introducing the series expansion of the Bessel function 

Xx 



F(x,s) 



dX B\x 2 (*- 2 ) 



Tl = 



_2(z~2) 
dX B\ 



(-iy 



oo 



(-iy 



n =o n\T {n+ 



Xx 



r in 



2 

z-2 



A 



2(z-2] 



2{z-2) 



-Xs 



-As 



we can identify the coefficients u n of the generating function 



u n {s) 



i-iy 



r n+-^a [2(*-2) 



l l 

z-2 2 JO 



dX B^X' ] 



-As 



(CIO) 



(C.ll) 



We note that, since s e [— 1; 0], the convergence of the integral requires that B\ decays 
sufficiently fast (at least exponentially) . In the following, we shall relate the B\ to the 
initial conditions u n (0). First, note that the initial time coefficient u n (0) is given by 

(-1)" 



U n (0) 



dX B x X r 



(C.12) 



Expanding now the exponential in the time-dependent coefficients (|C.11[) and 
identifying u n (0), it follows 

(-1)" 



u„{s) 



Tin 



2 

z-2 



[2(z-2)] 



i i 

z-2 2 JQ 



dX B x X 



n+- 



E 

.fc=0 



{-Xsl 
k\ 



(-1)' 



— fc!r (n+^J [2{z-2)] n+ -= 
[2(z-2)} kT ( n + k + lh) 



i _ i 

-2 2 JO 



dX B x A" 



E 

fc=0 



fe! 



u n+k (0)s k 



(C.13) 



Because of the asymptotic relation x b ^ a T(x+a)/T(x+b) = l+0(l/x) |U Eq. (6.1.47)], 
in the limit z — > 2 one recovers the result (|22j) or the linear chain. The final density 
of the n-string is obtained after performing the transformations (|16[) and (|C.2p and 
since Cq(0) = u n (0), one arrives at eq. (|40|) in the main text. 
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Appendix D. Number of stationary states in the two-species model 



Generalising the discussion of the single-species model [16j El 1321 G3], we briefly 
outline the analogous computation of the number of stationary states on a chain of 
L sites in the two-species pair-annihilation model of section 4. The counting is based 
on the observation that stationary states cannot contain neither AB nor BA pairs on 
nearest-neighbour sites. 

First, consider an open chain. Let Sl be the number of stationary states for a 

chain of L sites, and let Al := N ^ A \ L- 1 f ) denote the number of stationary states 
on a chain of L sites where the leftmost site is occupied by an A-particle. Similarly 
define B L =N (b \ L-1 f ) 



and O l := N 

small lattices, one has So := 1 and Si = 3, S-2 = 7, S3 = 17, £4 = 41 . 
contributions to Sl satisfy the recursions 



L-1 1 ) , such that S L = A L + B L + O l . For 

The three 



(0[ 



Or =N L-1 = N 



f 00 | L-2 f ) + N (®A \ L-2 f ) + N (®b\I^2 



= L _ 1 + A L _ 1 
= N ( A \ L-1 



B L _i — Sl-i 



B l 



N B\ L-1 



N(A0 
n(b0 



L-2 



L-2 



N [AA 
- N (bB 



L-2 



L-2 



L -1 + A L -1 

= Ol-! + B L -i 



Si 



(D.l) 
(D.2) 

(D.3) 
(D.4) 

(D.5) 

Second, consider a periodic chain. The number of stationary states is denoted by 

( A \ L-1 1 ^ denotes the number of stationary 



and collecting terms, we have the recursion 

Sl = 2Sl-i + Ol-i = 2Sl-i + Sl-2 
The solution is readily found, for all L > 



i(l 



L+l 



ill 



L + l 



1.207 (1 + V2 



S L = A L + B L + O l , where A L := N [A L-1 



states on a periodic chain of L sites where one of the particles is of species A and the 
trait indicates the periodic boundary conditions. Bl and Ol are defined similarly. 

For small lattices, one has Si = 3, S% = 7, S3 = 15, S4 = 35, Now, if one of the 

sites is empty, this opens the chain so that the number of stationary states is equal to 
the one of an open chain of L — 1 sites, hence Ol = Sl—i = Ol- Furthermore 

A L = N (i4|L-l| ) = n( A| L-2 |0lH + N ( A 



L-2 



A 



= N Ua\l^2 
= Al-i + N ( / 



N AA\ L-2 



= A 



L-1 



A 



L-1 



(D.6) 



since for the counting of the stationary states the pair AA acts in the same way as a 
single A-particle. Similarly, Bl = Bl-i +Bl—x- Adding these contributions, we find 

S L = S L -i + Al—i + B L -i + ~A~l—i + B L -i 

= Sl-i + 2 (Sl-i — Sl-2) 



1+ (1 + V2 



1 - V2 



1 + V2 



(D.7) 



where the expression in the last line, valid for L > 1, is obtained by using (|D.5l) . 

As found previously for the single-species model [TBJ [Ml [33], the number of 
stationary states is on an open chain about 20% larger than for a periodic chain. 
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